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ABSTRACT 

We present the first extensive and detailed theoretical scenario for the interpretation 
of Cepheid properties observed in the SDSS filters. Three sets of nonlinear connec- 
tive pulsation models, corresponding to the chemical compositions of Cepheids in the 
Milky Way, the Large Magellanic Cloud and the Small Magellanic Cloud respectively, 
are transformed into the SDSS bands by relying on updated model atmospheres. The 
resulting observables, namely the instability strip boundaries and the light curves, as 
well as the Period-Luminosity, the Wesenheit and the Period-Luminosity-Colour rela- 
tions, are discussed as a function of the metal content, for both the fundamental and 
the first overtone mode. The fundamental PL relations are found to deviate from linear 
relations when computed over the whole observed Cepheid period range, especially at 
the shorter wavelenghts, confirming previous findings in the Johnson- Cousins bands. 
The obtained slopes are found to be mildly steeper than the ones of the semiempirical 
and the empirical relations available in the literature and covering roughly the same 
period range, with the discrepancy ranging from ~ 13 % in u to ~ 3 % in z. 

Key words: Stars: oscillations — stars: variables: Cepheids 



1 INTRODUCTION 

Thanks to their characteristic Period-Luminosity (PL) 
and Period-Luminosity-Color (PLC) relations, Classical 
Cepheids are the most important stellar standard candles to 
estimate accurate distances within the Local Group, reach- 
ing ~ 30 Mpc from space observations. These Cepheid based 
distances provide an absolute calibration of the extragalac- 
tic distance scale and two HST projects have been devoted 
to this i ssue, obtaining an est i mate of the Hub ble constant 
(see e.g. iFreedman et alJfeOQll : ISaha et alj|200ll 1. Moreover, 
Cepheid properties are very important tools to investigate 
the physics and the evolutionary b ehaviour of intermediate 
mass stars (see e.g. Marconi 2009, and references therein). 
Even if the Cepheid PL and PLC relations are traditionally 
studied in the standard Johnson-Cousin BVRIJK magnitude 
system, the broadband niters e mployed in the Sloan D igi- 
tal Sky Survey (SDSS; see e.g. lAbazaiian et al.l I2003T ) are 
becoming more and more popular in the context of current 
and future projects. For example, observations with the VLT 
Survey Telescope (VST), the Canada- France-Hawaii Tele- 
scope (CFHT), the Pan-STARRSl survey (Kaiser 2004), the 
Dark Energy Survey2 (DES; Abbott et al. 2006), and the 
Large Synoptic Survey Telescope3 (LSST; Tyson 2002) are 
being done in the Sloan filters and from programs including 



a time-domain component, a large number of variable stars, 
including Cepheids, is expected. In particular, we are deeply 
involved in the INAF-GTO approved VST survey STEP 
(SMC in Time-Evolution of a P rototype interactin g late- 
type dwarf galaxy; PI: V. Ripepi. [Ripepi e~ al. 2006)). This 
project will include gsDSS and isdss time-series photome- 
try of the Bridge, connecting the two Magellanic Clouds, in 
order to study its stellar populations using the variable stars 
as tracers. Thanks to approved and submitted ESO proposal 
we plan to obtain 20-25 phase points in gsDSS and isdss on 
32 squared degrees to identify and study variables including 
both Classical Cepheids and RR Lyrae. Periods and mean 
gsDSS and isDSS magnitudes will be obtained from these 
data and accurate PL and PLC relations will be derived. 
In order to guarantee a reliable interpretation of these and 
other Cepheid data and to assure a comparison between the- 
ory and observations, in this paper we provide a theoretical 
scenario to interpret the properties of Classical Cepheids in 
the SDSS niters, for the characteristic chemical composition 
shown by these pulsators in the Milky Way and the Magel- 
lanic Clouds. 

The organization of the paper is the following. In Section [2] 
we present the nonlinear convective pulsation models de- 
veloped by our team, whereas in Section [3] the adopted 
model atmospheres to transform the theoretical quantities 
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and the bolometric light curves into the SDSS filters are dis- 
cussed. The results are shown in Section [4] and in Section [5] 
we present the comparison between model predictions and 
previous results available in the literature. The Conclusions 
close the paper. 



Table 1. Parameters of the computed pulsational models adopted 
in this paper. The last four columns are the temperatures in Kelvin 
of the first overtone blue edge (FOBE), the fundamental blue edge 
(FBE), the first overtone red edge (FORE) and the fundamental 
red edge (FRE) respectively. 



2 PULSATION MODELS 

The adopt ed nonlinear conve c tive p ulsati on models are 
taken from iBono et all (|l999al . l200ll . l2002f ) for the chem- 
ical compositions representative of Galactic and Magel- 
lanic Cepheids, namely Z =0.02, Y =0.28 (Milky Way), 
Z =0.008, Y =0.25 (LMC) and Z =0.004, Y =0.25 (SMC). 
For each selected stellar mass in the range from 3.5 M© to 
11 Mq a luminosity level based on an evolutionary Mass- 
Luminosity relati ons that neglects both core overshooting 
and mass loss (see lBono et al. l200ll . and references therein). 
Previous theoretical investigation has shown that the inclu- 
sion of mild overshooting (corresponding to an increase in 
the intrinsic luminosity of 0.25 dex, at fixed mass, according 
to Chiosi, Wood and Capitanio 1993) implies a decrease of 
about 0.2 mag in the distance moduls inferred from theoret- 
ical Wesenheit and PLC relations (see e.g. Caputo, Marconi, 
Musella 2000). 

The properties of adopted pulsation models are summarized 
in Table 1. For each combination of chemical composition, 
mass and luminosity, a wide range of effective temperatures 
has been explored to define the position of the instability 
strip boundaries for both the fundamental and the first over- 
tone mode. 

Thanks to their nonlinearity these models allow us to 
predict the full amplitude behaviour and the variations of 
all the relevant quantities along a pulsation cycle. More- 
over the inclusi on of a nonlocal time-dependent treatment of 
convection (seelStellingwerl 19821 ; IBono fc Stellingwerj[l993 ; 
IBono et aUll999al . for details) is the key ingredient to pre- 
dict the whole topology of the instability strip!]], as well as 
the detailed morphology of light, radius and radial velocity 
curves. The physical and numerical assumptions adopted in 
the computation of these mode ls are discussed in detail in 
IBono et ail (|l998l . Il999al . 120001 1 and are not repeated here. 
The modal stability has been investigated in the quoted pa- 
pers for both the fundamental and the first overtone modes, 
so that, moving from higher to lower effective temperatures, 
the following four edges are predicted for each selected chem- 
ical composition: the first overtone blue edge (FOBE), the 
fundamental blue edge (FBE), the first overtone red edge 
(FORE) and the fundamental red edge (FRE). Moreover, 
theo retical atlas of bolo metric light curves are available (see 
e.g. IBono et al l l2000bh for each combination of chemical 
composition and mass (luminosity), varying the model ef- 
fective temperature (period). 
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1 The existence of the red edge is related to the quenching effect 
of convection on pulsation. 



3 TRANSFORMATION INTO SDSS 
PHOTOMETRIC SYSTEM 

In order to obtain the pulsation observables of the in- 
vestigated Cepheid models into the SDSS bands, we 
have transformed both the individual static luminosi- 
ties and effective temperatures and the predicted bolo- 
metric light curves into the corresponding filters, us- 
ing model ATLAS9 non-overshoot ing model atmospheres 
|Castelli. Gratton fc Kurucj Il997al lbl). convolved with the 
SDSS transmission functions. 

In particular we used computed bolometric corrections and 
color indices, as tabulated for different effective tempera- 
tures (3500 K< Teff < 50000K), surface gravities (log g 
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from 0.0 to 0.5) and metallicities (\M/H]=0.5, 0.2, 0.0, 
0.5, 1.0, 1.5, 2.0 and 2.5, both for [a/Fe] = 0.0 and for 
[a/Fe] = 0.4 ), available on the homepage of Fiorella Castelli 
flhttp : / / wwwuser . oat . t s . astro . it /castelli /colors / sloaii.html [ ) . 
For each selected chemical composition, interpolation 
through these tables allowed us to obtain the static mag- 
nitudes, as well as the light curves and the resulting 
magnitude-averaged and intensity-averaged mean magni- 
tudes, in the u S dss, gsDSS, tsdss, isDSS, and z S dss filters 
(u,g,r,i,z in the following tables and figures) 



4 PREDICTED PULSATION OBSERVABLES 
IN THE SDSS FILTERS 

In this section we present the results of the transformation of 
pulsation model predictions into the SDSS filters, providing 
the first theoretical scenario to interpret observed Cepheid 
properties in these photometrical bands. 



4.1 The light curves and the Hertzprung 
progression 

The bolometric light curves have been transformed into the 
u, g, r, i, and z filters. In Figs. 1 we show some examples 
of the resulting light curves in the g filtefl for Z = 0.004 
(M = 7.15M ), Z = 0.008 [M = 7.15M Q ) and Z = 0.02 
(M — 6.5M ). The behaviours shown in these plots con- 
firm that the Hertzprung progression |Hertzsprung| 1 19261 : 
iLedoux fc Walravenl ll95a ). that is the relationship linking 
the phase of the bump and the pulsation period, occurs at 
P x 11.5 d, 11.0 d and 9 d for models at Z =0.004, 008 and 
0.02, respectively (see e.g. iBono et al.ll200ol ; iMarconi et al.l 
2005, for an extensive discussion of theoretical predictions 
for this phenomenon and the comparison with observations). 



4.2 The instability strip boundaries 

The predicted instability strip boundaries reported in the 
last four columns of Table 1 have been interpolated (we use 
a linear regression for the FOBE and FORE and a quadratic 
relation for the FBE and FRE) and transformed into the 
SDSS filters. The results and shown in Fig. [2] as a function 
of the adopted metallicity. 

The left panel of this figure shows the behaviour of the 
extreme edges (FOBE and FRE) of the pulsation region 
when moving from Z=0.004 to Z=0.008 and Z=0.02. The 
right panel only refers to the fundamental mode (FBE and 
FRE). 

These plots indicate that also in the g vs g-i plane the 
predicted instability strip gets redder as the metallicity in- 
creases. In particular the nonlinear trend and the metallicity 
effect, shown in the right panel for the fundamental mode 
boundaries, are known to directly reflect into th e behaviour 
of the predicted fundamental PL relations (see lBono et al.l 
Il999bl ). In other words, on the basis of the behaviours noted 
in this figure, the PL relation in the g filter is expected to 
be nonlinear and metallicity dependent. 



-3.5 
at 

-4 
-4.5 

-3.5 

DC 

-4 
-4.5 




!' 1 1 I 1 1 1 1 I 1 1 1 1 I 1 'J 
P = 12.4d (5400K)) 



I I I I | I I I I | I I I I | I I I 
:P = 10.6d (5400K) : 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

:P = 8.3d (5400K) 



!' 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 
P=11.6d (5500K) 



+-H- 


Z = 0,004 


4++ 


+H+ 


-H-H- 


-H4+- 




1 1 II jl 1 1 











P=11.2d (5450K) 




Z = 0.008 1 
I I I I I I I I I I I I I I I I I I 



P = 8.9d (5300K) 




!' 1 1 I 1 1 1 1 I 1 1 1 1 I 1 
P=11.3d (5550K) 



P=10.9d (5500K) : 



I I I I | I I I I | I I I I | I I I 
:P = 9.5d (5200K) 



Figure 1. Theoretical atlas of light curves transformed in the g 
band for Z=0.004 and M=7.15M (up), Z=0.008 and M=7.15M 
(middle) and Z=0.02 and M=6.5M (down). The arrow marks 
the position of the Hertzprung progression. 
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Figure 2. Predicted instability boundaries of the whole instabil- 
ity strip (left panel) and for the fundamental mode only (right 
panel) in the intensity-averaged g versus g — i plane, for the 
adopted model metal abundances. 



2 The light curves in the other filters and for other stellar masses 
are available upon request to the authors 
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4.3 Color-Color relations 

The transformed multifilter light curves can be reported 
in different colour colour planes. In Figs. [3] we show the 
theoretical fundamental (F) and First Overtone (FO) 
pulsators in the r — i versus g — r, i — z versus r — i and 
g — r versus g — i diagrams for the adopted metallicities. 
On the contrary of what happen in the g — r versus u — g 
plane, due to the significant sensitivity of the u band on 
metal abundance, in color-color plane reported in the figures 
the relations are very narrow and the effect of metallicity is 
smaller. In particular, this effect becomes negligible for the 
g — i and g — r combination of colors suggesting that the 
comparison between theory and observations in this plane 
could be used to evaluate colour excesses. 
The linear regression through the intensity-averaged mean 
magnitudes of synthetic F and FO pulsators (see subsection 
I4.4|l provides the following metal-dependent analytical 
colour colour relations: 



g - r = 0.12 + 2.34(r - i) (a = 0.009) 

i - z = -0.04 + 0.62(r - i) (a = 0.002) 

g - r = 0.04 + 0.70( 5 - i) (a = 0.003) 
for Z=0.004. 

g - r = 0.12 + 2.35(r - t) (a = 0.008) 

i-z = -0.04 + 0.63(r -i) (a = 0.002) 

g - r = 0.04 + 0.70(5 - i) = 0.002) 
for Z=0.008. 

g - r = 0.12 + 2.44(r - i) (a = 0.008) 

i-Z= -0.04 + 0.67(r -i) ((7 = 0.003) 

g-r = 0.04 + 0.71(5 - i) {a = 0.002) 
for Z=0.02. 

4.4 Period - Luminosity relations 

We used the s ynthetic distr i bution s of fundamental models 
computed by ICaputo et al.l |2000) to derive the PL rela- 
tions for different photometric bands listed in Table [2] for 
the three assumed chemical compositions. We fitted with 
a linear regression Mj= a+b LogP the predicted distribu- 
tion of F pulsators for different assumptions on the period 
range, namely logP ^ 2.0, logP ^ 1.5, logP ^ 1.0 and 
logP ^ 1.0 to determine the three corresponding slopes: 
b a ii, b^i.B, b s hort and bi ong . The investigation of the differ- 
ence between bshort and bi ong is relevant in the context of 
current studies in the lit erature claiming the existence of a 
break at 10 days (see e.g. ICaputo et al.ll2000l: [ Marconi et al.l 
l2005l : ISandaee et al.ll2004l : iNgeow et al.l 2008) As shown in 
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Figure 4. PL relations (logP 1.5, see Table 2) in different 
SDSS photometrical bands for fundamental (solid lines) and first 
overtone (dashed lines) models with Z=0.004.The same figures 
for Z=0.008 and Z=0.02, are given in the on-line version of the 
paper. 



Fig. [4] the PL relations get steeper by increasing the fil- 
ter wavelength, in agreement with the well-known empiri- 
cal and theoretical results in John s on-Cousins photom etric 
bands jMadore fc Freedmanl I"l99ll ; ICaputo et alJl200Ch . as 
well as with the s emiempirical PL relatio ns in the SLOAN 
bands derived by l|Ngeow fc Kanburll2007h . 
Moreover, in agreement with the results by Bono (2010), the 
slopes of the relations at LogP 1.0 are steeper than those 
for longer periods, indicating a clear break of the PL relation 
at 10 days, more evident for the Magellanic Cloud chemi- 
cal compositions, in agreement with empirical suggestions. 
However, the amplitude of this effect decreases as the filter 
wavelength increases. In the case of first overtone pulsators, 
that are intrinsically limited to a shorter period range, linear 
relations over the whole covered period range are derived. 

The results reported in Table [2] clearly show the ex- 
pected flattening of the various fundamental PL relations 
when the metallicity increases. Again, this trend decreases 
as the wavelength increases. 

The smaller, negligible passing from Z — 0.004 to Z = 
0.008, sensitivity to metallicity of first overtone PL relations 
is due to the much narrower instability strip that is also 
hotter than the one of fundamental pulsators, and in turn 
less affected by changes in the convective efficiency produced 
by metallicity variations. 



4.5 The Wesenheit relations 

As the PL relation reflects the instrinsic width of the insta- 
bility strip and its application implies the availability of a 
statistically significant sample, as well as the need to dis- 



© 0000 RAS, MNRAS 000, 000-000 



Predicted properties of Galactic and Magellanic Classical Cepheids in the SDSS filters 5 




0.0 0.5 1.0 1.5 
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Figure 3. Fundamental and first overtone (smaller points on the left side of each figure) models of different metallicities in colour - 
colour planes. A reddening vector for Ay = 0.5 is also shown in each panel. 



Table 2. Predicted slopes of fundamental and first overtone PL relations in the SDSS filters. Columns from second to fourth report the 
slopes in three different assumptions on the period range, namely logP ^ 2.0, logP ^ 1.0 and logP 1.0 and logP ^ 1.5. The last 
column reports the slopes of first overtone PL relations over the whole period range covered by models. 



Band ^F,all bp, short ^F,long bp.^l.S ^>FO,all 

Z=0.004 

u -2.503(0.016) -3.134(0.016) -0.897(0.088) -2.752(0.014) -2.963(0.015) 

g -2.893(0.012) -3.362(0.015) -1.789(0.066) -3.065(0.012) -3.136(0.016) 

r -3.098(0.009) -3.440(0.011) -2.301(0.050) -3.222(0.009) -3.291(0.012) 

i -3.177(0.008) -3.469(0.010) -2.494(0.043) -3.281(0.008) -3.363(0.010) 

z -3.231(0.007) -3.494(0.009) -2.620(0.038) -3.326(0.007) -3.410(0.009) 

Z=0.008 

u -2.274(0.016) -2.875(0.016) -0.849(0.088) -2.498(0.015) -2.751(0.013) 

g -2.707(0.012) -3.131(0.014) -1.766(0.065) -2.855(0.012) -2.919(0.014) 

r -2.949(0.099) -3.257(0.011) -2.265(0.049) -3.056(0.009) -3.080(0.010) 

i -3.044(0.008) -3.309(0.009) -2.459(0.042) -3.137(0.008) -3.155(0.008) 

z -3.109(0.007) -3.346(0.009) -2.588(0.038) -3.192(0.007) -3.206(0.007) 

Z=0.02 

u -1.156(0.0232) -1.833(0.038) -0.361(0.070) -1.503(0.025) -2.908(0.024) 

g -1.973(0.015) -2.307(0.031) -1.558(0.047) -2.151(0.018) -3.101(0.024) 

r -2.396(0.011) -2.632(0.023) -2.092(0.035) -2.527(0.013) -3.299(0.018) 

i -2.562(0.10) -2.763(0.020) -2.306(0.030) -2.672(0.010) -3.390(0.016) 

z -2.680(0.009) -2.852(0.017) -2.460(0.027) -2.774(0.007) -3.449(0.014) 
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Table 3. The Wesenheit relations for fundamental and first over- 
tone pulsators with different metallicities. 

W=a +b logP 



LogP 



Figure 5. Wesenheit relations in different SDSS photometrical 
bands for fundamental (solid line) and first overtone (dashed line) 
models with Z=0.004 The same figures for Z=0.008 and Z=0.02 
are given in the on-line version of the paper. 



entangle reddening and metallici ty effects on distance de- 
term inations (see discussions in ICaputo et al.l 1 19991 . |2000| . 
l200lf ). the reddening free Wesenheit relation is often pre- 
ferred. The latter relation relies on an assumed reddening 
law and includes a color term similar, but not identical, to 
the Period-Luminosity- Color color term (see below), thus 
partially reduci ng the effect of the finite width of the insta- 
bility strip (see ICaputo et all 2000). The predicted Wesen- 
heit functions in the SDSS filter s have been built by adopting 
the color coefficients given by iGiradi et"afl (120041 1 and are 
reported in Table 3 for both the fundamental and the first 
overtone mode. These results clearly indicate the decrease in 
the intrinsic dispersion when passing from the PL to the We- 
senheit relations, i n agreement with previous empirical and 
theoretical results ( Madore fc Freed man 199ll; [Caputo et all 



l2000l : iRipepi et alj|20l3 ). The obtained Wesenheit relations 
for Z=0.004 are shown in Fig. 5. These plots confirm the 
linearity and the small intrinsic dispersion of the Wesenheit 
relations that make them powerful tools to infer Galactic 
and extragalactic Cepheid distances. 

We notice that in agr eement with previou s results in the 
Johnson-Cousins bands (|Caputo et alj|2000h . the metallic- 
ity dependence of the Wesenheit relations is reversed when 
passing from optical to near-infrared colors with the smallest 
effect for the g, g — r combination. 



4.6 The Period-Luminosity-Color relations 

The Wesenheit relations are reddening free but only par- 
tially correcting the effect of the finite width of the instabil- 
ity strip with the inclusion of their color term. However, the 
true relation that holds for each individual Cepheid is the 
Period-Color-Luminosity (PLC) relation, directly descend- 
ing from the combination of the period-density, the Stephan- 
Boltzma nn and the evolutionary Mass-Lum inosity relations 
(see e.g. iBono et all Il999bl ; iMarconil 120091. and references 
therein). The disadvantage is that both individual redden- 
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ing and multiband photometry are required in order to apply 
PLC relations to infer Cepheid distances and this is the rea- 
son why these relations are rarely adopted in the literature. 
In any case the theoretical SDSS PLC relations based on 
the adopted pulsation and atmosphere models are reported 
in Table 4 and shown in Figs. 6, 7 and 8. We note that, dif- 
ferently from the PL and the Wesenheit relations, the PLC 
ones have been derived from the individual pulsating models 
and not from the simulated Cepheid distributions within the 
respective instability strips. This because of the mentioned 
intrinsic nature of these relations, expected to hold for each 
individual pulsator. Again we find a smaller metallicity ef- 
fect in the case of first overtone models. 



COMPARISON WITH PREVIOUS 
RELATIONS IN THE LITERATURE 



In a recent paper, iNgeow fc Kanburj (|2007h have derived 
semiempirical Cepheid PL relations in the SDSS magnitudes 
by combining observed LMC BVI mean magnitudes with 
theoretical bolometric corrections. In the same paper these 
authors construct empirical g,r PL relations, using the pub- 
licly available Johnson-Sloan photometric transformations, 
finding a good agreement with the semiempirical ones. 

In Fig. 9 we show the comparison between the PL slopes 
derived in this paper for fundamental models with Z = 0.008 
and logP < 2.0 (fi l led sy mbols) and the semiempirical ones 
jNgeow fc Kanburj |2007| . open symbols) for the it, g, r, i, z 
filters . In the case of g and r the empirical relations by 
the same authors (asterisks) are also shown. These compar- 
isons suggest that the slope of the PL relations obtained in 
this paper, using suitable pulsation and atmosphere models, 
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Figure 6. PLC relations for F (closed circles) and FO (open 
circles) models with Z=0.004. See also Table [4] 
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Figure 7. As for Fig. ©but for Z=0.004 
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Figure 8. As for Fig. ©but for Z=0.02 



Table 4. PLC Relations for fundamental and first overtone pul- 
sators with different metallicities. 
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for periods shorter than 100 d, are systematically steeper 
than the semiempirical and empirical ones presented by 
iNeeow fe Kanburl (|2007T l. with a deviation decreasing from 
~ 13 % in the u filter to ~ 3% in the z filter. 



6 CONCLUSIONS 

In this paper we have provided the first theoretical scenario 
to our knowledge for the interpretation of Classical Cepheid 
properties observed in the SDSS filters. Extensive and de- 
tailed sets of nonlinear convective pulsation models com- 
puted for the typical chemical compositions of Galactic and 
Magellanic pulsators have been transformed into the SDSS 
filters by using bolometric corrections and color-temperature 
transformations based on updated model atmospheres. The 
main results of this investigation can be summarized as fol- 
lows: 

• The transformed instability strip confirms the nonlin- 
earity of the fundamental boundaries and the dependence 
on the metal content already found in the Johnson- Cousins 
bands. This behaviour explains why PL relations are not 
linear at large periods. 

• The PL relations in the SDSS filters show evidence of 
the break at 10 days, in particular for the lower metal con- 
tents and the shorter wavelengths, as expected on the basis 
of previous empirical and theoretical studies. 

• The theoretical light curves in the SDSS filters present 
morphological characteristics similar to the ones in the 
Johnson-Cousins bands, with evidence of the Hertzprung- 
progression phenomenon at the expected period for each se- 
lected chemical composition. 

• The derived Wesenheit and PLC relations have, as ex- 
pected, a smaller intrinsic dispersion than PL relations and 
show a metallicity dependence that resembles the behaviour 
of their Johnson- Cousins counterparts. 

• The obtained relations for first overtone pulsators are 
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Figure 9. Comparison between the slopes of the theoretical PL 
relations de rived in this paper and t he semicmpirical and empir- 
ical ones bv lNeeow fc Kanburl d2007t) . 



all linear and much less affected by metallicity than the fun- 
damental mode ones. 

• A comparison between the relations obtained in this 
pa per and the semiempiric al and empirical ones obtained 
by iNgeow fe Kanburl (j2007T ) indicates a discrepancy in the 
slopes ranging from ~ 13 % in u to ~ 3% in z, with the 
theoretical relations being systematically steeper. 
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